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In Volborthite, spin-1/2 moments form a distorted Kagome lattice, of corner sharing isosceles 
triangles with exchange constants J on two bonds and J' on the third bond. We study the prop- 
erties of such spin systems, and show that despite the distortion, the lattice retains a great deal 
of frustration. Although sub-extensive, the classical ground state degeneracy remains very large, 
growing exponentially with the system perimeter. We consider degeneracy lifting by thermal and 
quantum fluctuations. To linear (spin wave) order, the degeneracy is found to stay intact. Two com- 
plementary approaches are therefore introduced, appropriate to low and high temperatures, which 
point to the same ordered pattern for J' > J. In the low temperature limit, an effective chirality 
Hamiltonian is derived from non-linear spin waves which predicts a transition on increasing J'/ J, 
from v3 x y/3 type order to a new ferrimagnetic striped chirality order with a doubled unit cell. 
This is confirmed by a large-n approximation on the O(n) model on this lattice. While the saddle 
point solution produces a line degeneracy, 0(1/ n) corrections select the non-trivial wavevector of 
the striped chirality state. The quantum limit of spin-1/2 on this lattice is studied via exact small 
system diagonalization and compare well with experimental results at intermediate temperatures. 
We suggest that the very low temperature spin frozen state seen in NMR experiments may be related 
to the disconnected nature of classical ground states on this lattice, which leads to a prediction for 
NMR line shapes. 



I. INTRODUCTION 

The study of frustrated magnetic insulators has wit- 
nessed a resurgence in recent times, with the discov- 
ery of a number of interesting materials with frustrated 
spin interactions. Amongst the most geometrically frus- 
trated lattices are the pyrochlore and the Kagome lat- 
tice, and perhaps the most interesting class of sys- 
tems are those that combine strong quantum fluctua- 
tions with frustration. Recently, spin-1/2 systems on 
the Kagome lattice have been identified, the mineral Vol- 
borthite Cu 3 V 2 07(OH)2 • 2H 2 Oi and Herbertsmithite^. 
In the former the equilateral Kagome triangles are dis- 
torted into isoceles triangles, rendering two of the near- 
est neighbor exchange constants different from the third. 
In the latter case, the Kagome lattice is believed to be 
structurally perfect, but with perhaps a small percent- 
age of impurity spins. Nevertheless, both systems display 
low temperature physics very different from their unfrus- 
trated counterparts, and do not show signs of ordering 
down to temperatures well below the exchange coupling 
strength. 

While a lot of theoretical effort has gone into charac- 
terizing the ideal frustrated lattices, distortions of the 
ideal structure, although common, have received less 
attention 3 -. In many frustrated magnets, lattice distor- 
tions occur spontaneously to relieve the frustration, lead- 
ing to a strong coupling between magnetic and struc- 
tural order parameters. Such 'multi-ferroic' couplings are 
highly prized from the technological viewpoint for conve- 
nient manipulation of magnetism 4 - and certain frustrated 



magnets are natural candidates^. This provides further 
motivation for studying the effect of distortions. From 
the theoretical viewpoint, the partial lifting of degeneracy 
from distortions can lead to a more tractable level of frus- 
tration, and allow for new theoretical approaches. Here, 
we consider the effect of lattice distortion on the Kagome 
lattice. The class of lattice distortions we focus on are 
motivated by the material Volborthite, whose structure 
consists of corner sharing isoceles triangles. Bonds along 
two directions then have exchange constant J while the 
bond along the third direction has exchange constant 
J' = aj. In Volborthite, it is not definitively known if 
a > 1 or a < 1, although a comparison of bond lengths 
seems to favor the former—. Hence we treat both kinds 
of anisotropy in this paper, with slightly more emphasis 
on the a > 1 case. 

We attack this problem first from the classical zero 
temperature limit. We show that for a wide range of dis- 
tortions, the large classical degeneracy of the Heisenberg 
model on the isotropic Kagome lattice is partially lifted, 
and the number of coplanar ground states now scales in 
a sub-extensive fashion, as the exponential of the lin- 
ear system size. An interesting comparison here is with 
the isotropic Kagome and pyrochlore lattices, where the 
extensively degenerate ground state can be specified in 
terms of local constraints reminiscent of the Gauss law 
of a lattice gauge theory^. In fact, that analogy has 
been carried further to describe new quantum phases of 
frustrated magnets corresponding to the coulomb phase 
of the lattice gauge theory*^. In contrast, the subexten- 
sive classical degeneracy of the distorted Kagome lattice 



2 



is naturally thought of as arising from constrains on one 
dimensional structures, and the 'soft-spin' dispersion on 
this lattice features a line degeneracy in the Brillouin 
zone, unlike the flat band of the Kagome lattice. Both 
these features are shared by pure ring exchange models on 
the square lattice as studied in Ref. IllL where a new spin 
liquid phase, the excitonic bose liquid, was discussed. 

In contrast to the isotropic Kagome system, the ground 
states of the distorted Kagome lattice are not connected 
by local moves, requiring moving an infinite number of 
spins to make transitions from one configuration to an- 
other. We suggest that this difference may be related 
to the experimental observation of spin freezing seen in 
NMR experiments at low temperatures in Volborthite 
(but not in the isotropic Kagome compound Herbert- 
smithite) . The classical ground state ensemble may then 
be expected to capture aspects of this glassy state, which 
we use to make experimental predictions. 

Next, we consider the question - if a spin system on 
this lattice develops long range magnetic order, what is 
the preferred spin pattern? The degeneracy is expected 
to be broken by fluctuation effects, and hence we ana- 
lyze the effect of quantum and thermally excited spin 
waves in the harmonic approximation. Remarkably, the 
spin waves are found to have a precisely flat dispersion, 
as in the ideal Kagome case, and do not distinguish be- 
tween the classical ground states at this level. To proceed 
we consider thermal fluctuations in the classical model 
with a > 1 in two complementary ways, approaching 
from high and low temperatures. These are found to 
be consistent with one another and point to a new fer- 
rimagnetic state, characterized by alternating chirality 
stripes, and a doubled unit cell, which we call the chiral- 
ity stripe state. The first calculation consists of combin- 
ing the low temperature non-linear spin wave expansion 
with the effective chirality Hamiltonian technique pio- 
neered by Henley 12 . While at the isotropic point our 
method picks out the -\/3 x \/3 state, consistent with 
many other studie o 12 ' 13 i 14 ' 15 i 16 , turning up the spatial 
anisotropy leads to a transition into a new state - the 
chirality stripe state. To attack the problem from the 
opposite, disordered limit, we consider generalization to 
the classical 0(n) model which is tractable in the limit 
n — > oo, and captures the fluctuating nature of the spins 
at high temperatures. At the saddle point level, the flat 
band degeneracy of the ideal Kagome case is shrunk down 
to a line degeneracy for a > 1 . Fully lifting the degen- 
eracy requires going to the next order in 1/n, which we 
accomplish by utilizing the high temperature expansion. 
The selected state has the same nontrivial wavevector as 
the chirality stripe state providing additional confirma- 
tion. In contrast, when a < 1, the large n saddle point 
itself picks out the q — wavevector. 

Finally, we study the problem in the quantum limit, via 
exact diagonalization studies on small (12-site) systems 
with spin- 1/2. Bulk properties such as specific heat and 
magnetic susceptibility at intermediate to high tempera- 
tures are found to be rather insensitive to the anisotropy 



and differences arise only below temperatures of about 
J/5, as seen in experiments^. On the other hand, the 
ground state of the small cluster is found to be a spin sin- 
glet and the spin gap decreases on increasing anisotropy. 

The structure of this paper is as follows. In Section HT1 
we discuss the classical ground states of the distorted 
Kagome model as well as the properties of the ground 
state ensemble, and possible connections to the NMR ex- 
periments on the low temperature state in Volborthite. 
Next, we address the question of which spin ordered 
pattern is favored by fluctuations on this lattice using 
two approaches, first by deriving an effective chirality 
Hamiltonian from non-linear spin waves in Section IIIII 
and next via a classical large-n 0(n) approach, in Sec- 
tion IIV1 which produce consistent results. Finally, the 
problem is treated in the extreme quantum limit via ex- 
act diagonalization of small systems in SectionfV] Details 
of calculations are relegated to three appendices. 



II. CLASSICAL GROUND STATES 

Consider the antiferromagnetic Heisenberg model on 
the distorted Kagome lattice (Fig. [I} with different cou- 
plings for bonds on the three principal directions, 



H 
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triangles 
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where S are quantum or classical spins, TJ J means 
Jab Jbc Jca, and A,B, C are indices for the three sub- 
lattices. 

If Jab, Jbc, Jca are all different, we call the lattice 
fully distorted Kagome lattice. In this paper however we 
consider mainly the distorted Kagome lattice in which 
Jab = Jca ^ Jbc- For simplicity we set J A b = Jca = 1 
and Jbc = a - The Hamiltonian simplifies to the follow- 
ing form, 



H= ^ ( s a • S B + aS B ■ S c + S c • S A ) 

triangles 

= f [(V^Sa + S B + S c ] 2 - constant 

triangles 



(i) 



There are two simple limits. In one limiting case, a — > 
0, the lattice becomes a decorated square lattice, with 
additional sites at the midpoints of square lattice edges. 
In the other, quasi-lD, limit a — > oo the lattice turns 
into decoupled antiferromagnetic chains and 'free' spins. 

From the lattice structure^, especially the Cu-0 bond 
lengths data, of Volborthite we expect that a > 1 in 
that material, although there is no direct experimental 
data available or quantitative first principles calculations 
available yet. Hence, the a < 1 case is also considered in 
some of the following theoretical treatments. 
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FIG. 1: (Color online) a L x L (L — 2) distorted Kagome 
lattice. ei,e2,e3 are three lattice translation vectors. The 
exchange constant J for bonds along the e2 and e3 direc- 
tions are equal, but different from J' = aj, the exchange 
constant for bonds along the ei direction. For Volborthite, 
it is believed that a > 1. T, M, K are high symmetry points 
in the Brillouin zone (BZ). The proposed spin-ordered state 
has alternating positive/negative chirality stripes and Fourier 
components at T- and M-points in BZ. Dashed rectangle is 
the reduced BZ for the doubled magnetic unit cell. 



The first step of studying the classical ground states on 
the lattice is to solve the classical ground states of a single 
triangle. Setting the 'cluster spin' in Eqn. JT]) to zero we 
can solve the angle between A-site spin and B(C)-site 
spin, denoted as 9o = arccos(— l/2a) (see Fig. [2|. Since 
a =/= 1, this angle will be in general incommensurate to 
2ir. We ignore the accidental commensurate cases in the 
following discussion since they form a measure-zero set 
of a. Then the (3-state) Potts model description for the 
coplanar ground states of the isotropic Kagome case does 
not work for the distorted Kagome lattice. 

A special case is a < 1/2. In this case there is no way 
that the 'cluster spin' can be zero and the classical ground 
state is a collinear state with A-site spin anti-parallel to 
B(C)-site spin (Fig. [2]). Thus, for a < 1/2, the classi- 
cal ground state is collinear and there is no degeneracy 
except a global spin rotation. Notice that the lattice be- 
comes bipartite (not frustrated) in the limit a = 0. This 
classical consideration shows that the frustration of BC 
bonds is ineffective for nonzero a < 1/2. Later we will 
see from exact diagonalization study that this naive clas- 
sical picture survives in quantum regime. The classical 
collinear state has a macroscopic net moment and is a 
'ferrimagnetic' state. 

For a > 1/2 case we expect that coplanar classical 
ground states are favored by thermal or quantum fluctu- 




FIG. 2: (Color online) ground states of a single triangle 
(8o — arccos( — l/2a)), definition of chirality variables, and 
constraint on the six chirality variables for the distorted 
Kagome model, on a single hexagon of the honeycomb chi- 
rality lattice. ABCDEF are labels of six spin sites used to 
calculate the effective chirality interactions in a later section. 

ations, and there will be zero energy band(s) for the 0(n) 
model with n > 3, because the Hamiltonian in Eqn. ([1]) 
can be written as a sum of squares of 'cluster spins—. 
Then it is convenient to utilize the chirality variables used 
in the isotropic Kagome modei^. The chirality variables 
are Ising variables living at the centers of triangles, thus 
forming a honeycomb lattice. The positive or negative 
chirality variable represents the cluster of three spins on 
a triangle rotating counter-clockwise or clockwise when 
one goes from A- to B- then to C-site, or Sa X (Sb — Sc) 
pointing toward the +z or —z direction, assuming all 
spins lie in the x — y plane (Fig. [2]) . 

It should be emphasized that the chirality variables 
are not independent. They determine how spin rotates 
(counter-clockwise or clockwise) when one walks along 
a bond, but after walking along a closed loop on the 
lattice the spin should go back to the initial direction. 
We need only to consider length-six hexagonal loops on 
the (distorted) Kagome lattice. Each one of these loops 
will impose a constraint on the six chirality variables r\ 
in the corresponding hexagon in the honeycomb chirality 
lattice (Fig. HZ), 

■ni(29 )-7] 2 e -7] 3 e Q + T] i (2e )-T] 5 e -r] 6 9 = mod 2tt 

(2) 

For the isotropic Kagome antiferromagnet, #0 = 27r/3 
and the constraint simplifies to J2i=i % = ^ or 0. There 
are 22 allowed patterns on a single hexagon out of 2 6 = 64 
combinations. For the distorted Kagome model, 9q is in- 
commensurate to 27r and the constraint is more restric- 
tive: £)i=i r] t = ±6 or, £\ =1 rji = and rji + 774 = 0. 
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The last equation is the new constraint compared to the 
isotropic Kagome lattice. Note, this constraint holds for 
all a 7^ 1, so long as a coplanar ground state is favored, 
i.e. a > 1/2. There are only 14 allowed patterns on a 
single hexagon. For the fully distorted Kagome lattice, 
the constraint is even more restrictive: Sj=i Vi ~ ^6 or > 
T)i + Vi = V2 + r l5 = V3 + % = Q- There are only 10 
allowed patterns on a single hexagon. 



A. Properties of Coplanar Ground States on 
Distorted Kagome Lattice 

The degree of degeneracy for these models on a lattice 
is a much more subtle problem. From Baxter's solution^ 
we know that the degeneracy of the coplanar ground 
states of the isotropic Kagome antiferromagnet (or 3- 
state Potts antiferromagnet) is extensive, exp(0.379iV), 
where N is the number of Kagome unit cells. 

By counting the allowed chirality patterns for the dis- 
torted Kagome model with the L x L open boundary 
geometry in Fig. Q] up to L — 9, we conclude that the 
degeneracy is 'sub-extensive', about exp(2.2L). Table U 
lists the exact enumeration result. 

Appendix [A] derives the asymptotic formula of the de- 
generacy by transfer matrix method for a slightly dif- 
ferent geometry with periodic boundary condition. The 
'sub-extensive' behavior is proved by rigorous upper and 
lower bounds and the asymptotic formula. 

For fully distorted Kagome model the degeneracy is 
also 'sub-extensive', about exp(1.4L) for the geometry in 

Fig.ru 

One should be aware that the constant in the exponent 
depends on geometry and boundary conditions. Notice 
that Appendix lAl uses another geometry so that the re- 
sult is not exactly the same as the enumeration results, 
although they both show 'sub-extensive' behavior. 

Another issue about classical degeneracy is the ex- 
istence of the so-called 'weather-vane' modes. In the 
isotropic Kagome model those local zero-energy modes 
were argued to favor the V3 x y3 stated. However one 
can easily prove that in the distorted Kagome 0(3) model 
there is no local 'weather-vane' modes. This is because 
the cluster of spins of a 'weather-vane' mode must be 
bounded by spins pointing to the same direction. Those 
boundary spins inevitably involves all three sublattices if 
the cluster is finite. But an A-site spin can never be in 
the same direction as a B-site spin if 9$ = arccos(— l/2a) 
is incommensurate to 2n. 

There could still be non-local 'weather-vane' modes in- 
volving an infinite number of spins in the thermodynamic 
limit. But the number of these modes do not scale as 
the area of the system. In this respect, the ground state 
manifold of the distorted Kagome model is much less con- 
nected than that of the isotropic Kagome model. Thus 
glassy behavior is more likely to happen in the distorted 
model. 

Huse and Rutenberg studied the ground state ensem- 



ble of the isotropic Kagome antiferromagnet^ by field 
theoretical and Monte Carlo methods, and found that 
the spin-spin correlation has the V3 x V3 state signature 
but with power-law decay. 

We study the classical ground state ensemble of the 
distorted model by measuring the ensemble averaged 
spin-spin correlation. Lacking a good Monte Carlo al- 
gorithm we use the exact enumeration result for L x L 
lattice with open boundary up to I = 9. Because of 
the small size and possible boundary effects we have not 
been able to extract the scaling form of the correlation 
functions. However the result is qualitatively different 
from those of the isotropic Kagome antiferromagnet. For 
A-sublattice the correlation has a large q = (r-point) 
component. For B(C)-sublattice the correlation has a 
large Fourier component at the M-point, the mid-point 
of the BZ top(bottom) edge. 

Based on these hints we propose an ordering pattern 
as in Fig. [1] It has horizontal alternating stripes of posi- 
tive(negative) chiralities. We will later call it the chirality 
stripe state. This pattern doubles the magnetic unit cell 
in the vertical direction, thus reduces the BZ, and the M- 
point is actually equivalent to the r-point for the reduced 
BZ(Fig.[T|). 

To further confirm this we measured the mean-square 
of three Fourier modes of the chirality variables (m 2 ): 
(i) the uniform pattern, corresponding to the q = (T- 
point in BZ) spin configuration, with rnr = y^V, (ii) 
the staggered pattern, corresponding to the V3 x \/3 
spin configuration of the isotropic case or K-point in 
BZ, with tt«k = X) where the two sublattices in 
the honeycomb chirality lattice have opposite ± sign; 
and (iii) the chirality stripe pattern, corresponding to 
our proposed spin configuration (M-point in BZ), with 
m M — ±?yexp(ikM • R) where the ± signs are the same 
as the staggered pattern, R is the position of the honey- 
comb unit cell, kM is the wavevector of M-point (Fig.[l|. 

Results are summarized in Table HI For the isotropic 
Kagome model, the staggered pattern mode has the 
largest mean-square value, while for the distorted 
Kagome model the chirality stripe pattern has the 
largest mean-square value, which is consistent with the 
ensemble-averaged spin-spin correlation result. Also 
from the scaling of the mean-squares with system size 
we conclude that there is no long-range-order for chiral- 
ity variables at these Fourier modes. 

B. Comparison with 51 V NMR in Volborthite 

We have already noted that in contrast to classi- 
cal ground states on the isotropic Kagome lattice, all 
ground states on the distorted Kagome lattice are dis- 
connected from one another, and require moving an in- 
finite number of spins. Within a semi-classical view- 
point, large kinetic barriers separating the distorted 
Kagome ground states might lead to freezing at low 
temperatures. Interestingly, low temperature NMR 
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TABLE I: Exact enumeration results for L x L open boundary 
chirality lattice in the geometry of Fig. [T] The number of clas- 
sical ground states -/Vgs for isotropic and distorted Kagome 
lattices are shown. The tendency to order in different patterns 
[q = (F), \/3 x \/3 (K) and stripe (M) patterns] are com- 
pared by evaluating mean-square values of relevant chirality 
combinations [(mf), {m%), (m^) respectively]. 





Kagome 




distorted 




L 


AT < m r> 
NGS <^> 




A^GS 




( m M> 


1 


22 0.50 


1.00 


14 


0.64 


1.00 


2 


952 0.32 


0.90 


168 


0.62 


1.03 


3 


84,048 0.22 


0.92 


1,864 


0.61 


1.25 


4 


15,409,216 0.17 


0.84 


19,724 


0.61 


1.25 


5 






201,584 


0.61 


1.31 


6 






2,008,276 


0.61 


1.35 


7 






19,596,536 


0.61 


1.45 


8 






188,078,644 


0.60 


1.41 


9 






1,779,795,056 


0.60 


1.48 



experiments** on Volborthite indicate spin freezing be- 
low 1.5K(~ J/60)^, but no such freezing is observed in 
the isotropic Herbertsmithit o 21 i 22 ' 23 i 24 . It is tempting to 
attribute this difference in behavior to the difference in 
connectivity of classical ground states in the two cases. 
The vanadium atoms occupy the hexagon centers of the 
Kagome lattice, and are hence coupled to six spin-1/2 Cu 
moments. Experimentally, on cooling through the glass 
transition temperature there is a rapid rise of 1 /T% , and 
at lower temperatures two distinct local environments 
for the 51 V sites appear, a higher static field environ- 
ment (rectangular lineshape) estimated to involve 20% of 
spins, and a lower field environment (gaussian lineshape) 
for the remainder. We assume that the glassy state lo- 
cally resembles one of the classical ground states, and 
that they occur with equal probability. Then a volume 
average of a local quantity in the glassy state corresponds 
to an ensemble average over classical ground states. Of 
relevance to the NMR experiments here is the distribu- 
tion of exchange fields at the 51 V site, arising from spin 
configurations on the hexagons. For the nearly isotropic 
case awl, three different field values (H) are possi- 
ble, H sa 3H CU , H sa \f~ZHcu and H ss 0, where Hc u 
is the field from a single spin. For example, the first 
corresponds to a local a/3 x v3 pattern with staggered 
chirality. We need to calculate the probability to find 
these different fields. 

The authors of Ref. [lj^ put forward the interesting sug- 
gestion that the high field component seen in NMR cor- 
responds to local a/3 x -\/3 pattern. Their arguments 
though rested on properties of the isotropic Kagome 
model. Here, we analytically evaluate the probability dis- 
tribution of different field configurations for the distorted 
Kagome lattice using the transfer matrix method (details 
in Appendix A). The probability of obtaining the 3H CU 
exchange field is found to vanish in the thermodynamics 



limit, while that of the ^/3Hc u is 25% and of the approx- 
imately zero field configuration is 75%. This is roughly 
consistent with the experimental observation, but implies 
a revised value for the local moment that was obtained 
in Ref. 19 which assumed a local field of 3Hc u - Hence we 
anticipate a copper moment per site of 0.4 x -s/3 = 0.7 of 
the full moment, for small anisotropy. If the anisotropy 
is significant, the local field also changes, with the pre- 
vious y/3Hcu — > \J (5a — 2) jo?Hc u and the zero field 
values now being |2 - 2a- 1 \H Cu (with 50% probability) 
and a _2 |a — l\H Cu (with 25% probability). This sug- 
gests an upper bound for the anisotropy by requiring the 
local moment be less than unity, which gives a < 1.6. 

III. EFFECT OF FLUCTUATIONS ABOUT THE 
CLASSICAL GROUND STATES 

It is well-known that thermal or quantum fluctua- 
tion can lift the classical ground state degeneracy- 2 ^. In 
the isotropic Kagome model these kinds of 'order-by- 
disorder' studies suggest that the Kagome antiferromag- 
net would select the a/3 x a/3 ground state^l 3 -^^^, 
namely the staggered chirality pattern. 

We study the 'order- by-disorder' effect in the distorted 
Kagome model (a > 1) by quantum and classical 'spin 
wave' theory. It is found that at quadratic order the fluc- 
tuations (quantum or classical) cannot distinguish dif- 
ferent coplanar classical ground states. One has to go 
beyond quadratic order of fluctuation to find 'order-by- 
disorder' phenomenon. 

A. Linear Spin Wave Theory 

A classical coplanar ground state can be described by 
angles 6j of classical spins with respect to a reference 
direction in spin space. Define a local spin axis for every 
site such that the S z axis is perpendicular to the common 
plane of all classical spins, and the S x axis is along the 
classical spin direction. 

The Hamiltonian becomes 

H=J2 MSfS^ + cosie^iSfS^ + SfS]) 

<H> (3) 
-sm(%)(Sf SV-SfSJ)] 

where #y = Bi — 9j is the angle between classical spins 
on sites i and j, and the chiralities determine the sign of 
these angle differences. 

For quantum spin-S* spins we can use the Holstein- 
Primakoff bosons to describe the fluctuations 

S? = Sf + iSf = V25 - m ■ h 
Sr = SV-iS? = bly/2S- ni 
where rii = b\b i is the boson number operator. 
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Expanding in powers of 1/S, the Hamiltonian becomes 

H = E GS + S' 3/2 H 1 +S-H 2 + S 1/2 H 3 + H 4 + ... 

where Eqs is the classical ground state energy, H n con- 
tains n-th order boson creation(annihilation) operator 
polynomials. In fact Hi identically vanishes. H2 gives 
the quadratic (or so-called 'linear') spin wave theory. 

2 = £ H 2>i j 
<ij> 

H 3tij = - Jij cos(%)K + nj - (l/2)(fot + &.)(&} + b s )] 
-(l/2)J«(6j- 60(6} -6,.) 

(4) 

Notice that H2 only depends on cos(#y), then it is iden- 
tical for all classical ground state configurations can 
differ only by a sign between different classical ground 
states). Therefore spin wave expansion at the quadratic 
level cannot lift the classical degeneracy. 

Dispersion of the quadratic spin wave is presented in 
Appendix [Bl One interesting result is that although the 
dispersion becomes much more complicated than that of 
the isotropic Kagome model, the zero-energy fiat band 
still exists. Another strange feature is that as long as 
q^I, a > 1/2, the 'spin wave velocity' vanishes in the 
direction perpendicular to the BC-bonds. 



B. Classical 'Spin- Wave' Expansion and Effective 
Chirality Hamiltonian 

To lift the classical degeneracy we need to consider the 
'non-linear' spin wave theory, especially the cubic order 
terms H 3 , because they are the lowest order terms dis- 
tinct for different classical ground state configurations. 
Following Henley and Chani 2 - we can in principle derive 
the effective interactions between chirality variables. In 
the remaining part of this section we use a different for- 
malism by combining Henley's idea and the classical low 
temperature 'spin wave' expansio n 27 ' 28 ! 29 . 

We consider classical 0(3) spins on the distorted 
Kagome lattice. To simplify the notations we set the 
spin length S to unity. We define local spin axis as in 
previous subsection, S z axis perpendicular to all spins, 
S x axis along the classical spin. We can still use the ex- 
pression Eqn. ([3]) for the Hamiltonian. For classical spin 
it is convenient to parametrize the fluctuation by 



and the in-plane e y and out-of-plane e z fluctuations are 
supposed to be small at low temperatures. 

The most important contributions to the partition 
function comes from fluctuations around classical ground 



states. 

Z = Z o 1 J 2?Sexp(-/?tf)]p[(S 4 ) 2 -1] 

% 

ex J2 [veyVe z exp(-pH)l[(l/S?) 

classical GS i 

where <5[(Si) 2 — 1] is the Dirac-(S function used to ensure 
unit spin length, the product Yl{l/Sf) is the Jacobian 
of changing variables from 0(3) spin to e v and e z . Z = 
(2tt) 3N is chosen in such a way that Z — > 1 as — ► (N 
is the number of unit cells). 

Absorb the Jacobian into the exponential and expand 
5"" in terms of e v and e z , then the exponent becomes 

-0H = -0(H% + H z +H 3 + H 4 - (1/2)T^& + . . . ) 

i 

where — (1/2)T^ [(ef) 2 + ( e f) 2 ] comes from the Jaco- 
bian, and to simplify the notation we define £j = {c() 2 + 
(ef) 2 . Then 

HI = ]T J tj cos(%){e^ - (l/2)[( e f) 2 + (ef) 2 ]} 

<tj> 

HI = £ - (1/2) Jij cos(%)[(6f) 2 + (e z f] 

(5) 

ff 3 = (l/2)££^sin(%)e^ 

i 3 

H 4 = (1/8) £ J l3 cos(%)(^ 2 + + i 2 ) 

<ij> 

Again the quadratic terms are identical for all classical 
ground states. 

We can rescale e v and e z by y/0 to absorb (3 into H\ 
and iff- Define e y = yf]3e v ,e z — \/~j3e z , then the expo- 
nent becomes 

-0H= -H$ - H z 2 - VTH 3 - TH 4 - 0(T 2 ) 

where H%' z , H 3 are obtained by replacing e y , e z by e y , e z 
in the formulas of H^.H^H^, respectively. H4 com- 
bines the original quartic order term H4 and the lowest 
order term from the Jacobian, and we have set the Boltz- 
mann constant — 1. Since higher-than-quadratic or- 
der terms are controlled by temperature, we can do a 
controlled perturbative expansion in powers of the small 
parameter T. 

As the first approximation we may keep only H\ and 
-fff for very low T. Solution of the quadratic theory 
is presented in Appendix [Cj The out-of-plane fluctua- 
tion e z has a flat zero-energy band, which is consistent 
with Moessner and Chalker's mode-counting argument^. 
The in-plane fluctuation has the 'Goldstone' mode at 
wavevector k = 0. But since this is a classical theory, 
the dispersion around the 'Goldstone' mode is quadratic. 
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C. Effective Chirality Hamiltonian 

Now we can formally write down an expansion for small 
T. Define Z = J exp(-H% - H z )Vi y Vi z . Remember 
that Z is the same for all classical ground states we are 
perturbing. The free energy / per unit cell for fluctua- 
tions around one classical ground state is 

/ = (l/N)E GS — 3T In T — {l/N)T\nZ 

- (1/2)T 2 ((H 3 ) 2 /N) Q + T 2 (H i /N) Q + 0(T 3 ) (6) 

where N is the number of unit cells, (A) means the 
expectation value in the quadratic theory, i.e. (A) = 
Zo 1 J A ■ exp(-iTf - H%)Vl y Vl z . Since Z and Hi are 
identical for all classical ground states, difference at T 2 
order comes from the ((H3) 2 / N) term. Remember that 
each term in H3 contains a sin(^y), the sign of which is 
determined by the chirality of the triangle containing the 
bond < ij >. Therefore ((H 3 ) 2 /N} will generate effec- 
tive chirality-chirality interactions Jijijir/j for each pair 
of chirality variables rji and rjj . Details about calculating 
the chirality interactions are presented in Appendix [Cl 

There are two technical obstacles for this 'order-by- 
disorder' analysis: (i) The flat zero-energy band will 
make the two-e z correlation function diverge; (ii) The 
'Goldstonc' mode will make the two-e y correlation func- 
tion diverge. Both (i) and (ii) will make ((H 3 ) 2 /N} 
divergent. 

To proceed we add a term J z J2i {Si) 2 m the Hamil- 
tonian. This can be thought as a single-ion anisotropy 
term disfavoring out-of-plane fluctuation. The flat zero- 
energy band will be shifted to a positive value and no 
longer produce divergence. We also need to cure the di- 
vergence from the in-plane 'Goldstone' mode. But no 
natural interaction can do this job. Therefore we add 
an artificial mass term J v (ef ) 2 to the Hamiltonian, 
which gives the 'Goldstone' mode a small gap, or can be 
thought as an infrared cutoff. Eventually we would like 
to take the limit J z , J y -> 0. 

To check consistency we first calculated the effec- 
tive chirality interactions for a = 1 Kagome model. 
The interactions are antiferromagnetic and seems to be 
short-ranged (see Table |H] in Appendix [C]). Because 
the nearest-neighbor chirality antiferromagnetic coupling 
dominates, the staggered chirality pattern (namely the 
spin configuration) is selected, which is consis- 
tent with all previous 'order-by-disorder' studies for the 
isotropic Kagome model. This selection is independent 
of J z and J y for the range of parameters we studied. 

The a > 1 case is more delicate. It seems that the 
chirality interaction is not short-ranged (see Table |TT] in 
Appendix [C]), and the selection of chirality pattern is 
more sensitive to J z and J y . We have calculated chiral- 
ity interactions up to sixth-neighbor, with J v = 0.01 as 
the smallest value we can use, and for various J z and 
a. A rough picture (Fig. \5§ is that for a close to unity 
or small J z the staggered chirality pattern (analogue of 
the V3 x spin configuration of the isotropic Kagome 



J z J Y =0.01 




a 



1.1 1.2 1.3 



FIG. 3: Phase diagram obtained from classical spin-wave 
'order-by-disorder' analysis with J y — 0.01 (artificial gap for 
in-plane fluctuation). J z is the single- ion anisotropy. For 
large a or J z , the chirality stripe pattern is selected. Dots 
are calculated boundary points and the line is a guide for the 
eye. Weak interlayer couplings are assumed, to stabilize true 
long-range order. 



model) is still favored, but in the other part of the pa- 
rameter space our proposed chirality stripe state is se- 
lected. One should be aware that this picture may still 
depend on the unphysical parameter J y , and including 
further neighbor chirality interactions may also modify 
the phase boundary 

We have also tried to use this approach for a < 1 
case. However the selection of ground state is much more 
sensitive to the unphysical parameters and the number 
of chirality-chirality couplings we include. We decide to 
leave this part for more detailed studies in the future. 



IV. LARGE-N APPROXIMATION 

Another way to study (anti)ferromagnet is to general- 
ize classical 0(3) spin to O(n) spin. At n — > 00 limit the 
theory can be solved exactly by saddle point approxima- 
tion. One can also calculate 1/n corrections systemati- 
cally The saddle point approximation is supposed to be 
good for high-temperature disordered phase. As temper- 
ature decreases one can usually decide at which wavevec- 
tor the long-range-order is developed, by looking at the 
position of the lowest 'excitation' energy, or the lowest 
eigenvalue(s) of the inverse of the spin correlation func- 
tion matrix. 

For the isotropic Kagome model the lowest 'excitation' 
is wavevector independent in the saddle point solution. 
For distorted Kagome (a > 1) model the lowest 'exci- 
tation' is degenerate on a line in momentum space. We 
have to include 1/n correction to determine the possible 
ordering wavevector uniquely. However, for a < 1 the 
q = wave-vector is selected at the saddle point level. 
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A. Saddle Point Solution and Line Degeneracy for 
Distorted Kagome Lattice 

The model we use is the 0(n) spin antiferromagnet on 
the distorted Kagome lattice. 

H = £ ^(515g + a-Sg5g + 5g51) (7) 

triangles a 

with constraints Y^l=i (^i) 2 = 1- We rescale all spins 
and /? by Sf = ^/nSf, (3 = (3/n. The partition function 
becomes 



id/ J 



Z = Z o' J ( II d ^ I exp(-^F) J] <5 

\i, a J i 



-£(sf) 5 



where Z = [?W 2 7r™/ 2 /r(n/2)] 3Ar such that Z -► 1 as 
/3 — > 0, iV is the number of unit cells, is the Hamilto- 
nian if with S directly replaced by S. In the remainder 
of this section we will write S, H and (3 as S, H and (3, 
respectively. We will write VS instead of JT^ a dS? ■ 
Using the fact that 



S(x) 



— exp + {li )a 



where Xi is a real dummy variable, and is an arbitrary 
real parameter to be determined later by the saddle point 
condition, we can rewrite the partition function as 

Z =Zq 1 J VSVX exp | -f3H + ^ A* ™ ~ £ ( S if j 



Zq L / PSPAexp 



(8) 



in which Mij = (tA; + + f3Jij/2 is a symmetric 

matrix, A = iX + fi, and PA = J\ i [dXi/{2w)]. Integration 
over Sf gives 



Z = Zo\ 3nN / 2 J VX det(M)- n / 2 exp[J]n(iA J 



Now the saddle point condition is 



d 
d\ii 



lndet(M)=2, Vi 



Let us assume the saddle point solution has all lattice 
symmetries, e.g. translational invariance. Then ^ de- 
pends only on which sublattice the site i belongs to. Fur- 
thermore, because the B- and C-sublattices are equiva- 
lent, we have /ib = Mo 

Assuming translationally invariant /Kj, the matrix 
can be block-diagonalized by Fourier transformation. 



\KM 
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a=1.5] 








M 


r 


k: m: 




FIG. 4: (Color online) dispersion u of the O(n) model in 
the saddle point approximation along certain high symmetry 
directions (shown in the first panel), for three different a, 
at P = 0.2. The lowest eigenvalue(s) are shown with bold 
lines(dot). 



det(M) = n k det[Af(k)] where M(k) is a 3 x 3 matrix 

' fi A /3cos(V2) /?cos(/c 2 /2) 
M(k)= /3cos(fc 3 /2) (is a/3 cos(fci/2) 
y/3cos(fc 2 /2) a/3cos(fci/2) /ic 

with fci = k • (/c3 = — fci — k,2). The saddle point 
condition becomes 



(1/A0 



d 



lndet[M(k)] = 2, A = A,B,C 



and in the thermodynamic limit N — > oo the sum 
becomes a integral over Brillouin zone, (l/^OX^k — > 
f f dfc l( ifc 2 /(2,) 2 

This saddle point equation cannot be solved analyti- 
cally. But when (3 is small, we can expand it in terms 
of (3 and obtain a high-temperature series for \ix- The 
result is 

2/i A = I + 4[3 2 - 4a(3 3 + . . . (9a) 
2MB.C = 1 + 2(a 2 + l)/? 2 - 4a/3 3 + . . . (9b) 

We notice that /ib,c > Ma for a > 1, which leads to a 
degenerate line of lowest excitation in the saddle point 
approximation. This high-temperature (small (3) series 
can be extended to intermediate temperature {(3) by Pade 
approximation. 

After solving fix we can solve the 'dispersion', or the 
eigenvalues of M(k). Dispersion along certain high sym- 
metry directions are shown in Fig. [4] Note, for a < 1 
the lowest eigenvalue is uniquely determined at k = 0. 
However, for a > 1 the lowest eigenvalue is degenerate 
on the ki — line, or the vertical T — M line in the BZ. 
Finally, for a = 1 the lowest eigenvalue is degenerate over 
the entire BZ. 

To decide the ordering wave-vector uniquely we must 
consider 1/n correction for a > 1 cases. Before 
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FIG. 5: (Color online) predicted elastic neutron scattering in- 
tensity for distorted Kagome model, obtained from the saddle 
point solutions at /3 = 0.2 for four different a. Hexagons are 
BZ borders. Darker region have lower intensities. The a > 1 
case shows quasi-lD feature. 
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presenting that in the next subsection, we show the 
calculated elastic neutron scattering intensity (Fig. O 
[J2x y (M~ 1 )xy] 2 of the saddle point solutions for four 
different a with relatively high temperature (3 = 0.2 
(summation is over X,Y = A, B, C). We emphasize that 
the maximum appearing in the elastic neutron scatter- 
ing intensity does not directly correspond to the possible 
long-range-order wavevector. 

B. Lifting the Line Degeneracy of a > 1: 1/n 
correction 

To lift the degeneracy of the lowest 'excitations' of the 
saddle point approximation, we have to include fluctua- 
tions around the saddle point. 

We have three Ax.k fields and 3n S x k fields in the 
action, where X is the sublattice index, a is the O(n) 
index of spin. The Green's function of the spins with the 
same O(n) indices is a 3 x 3 matrix. Under the saddle 
point approximation it is G^ ' abxy (k) = [A/(k)]^-y<5 a b 
where X, Y = A, B, C for three sublattices, a, b are O(n) 
indices. We need the correction to this Green's function 
by the fluctuations of A around zero. From Eqn. ([5]) we 
see that there is a three-leg vertex between A and S a , of 
the form -iXi(Sf) 2 . 

The Feynman rules and Dyson equations are summa- 
rized in Fig. [HI Notice that the three-leg vertex preserves 
sublattice index for all fields and also O(n) index for the 
spins. There is no free propagator for A fields in the orig- 
inal theory. To make the perturbative expansion well 
defined we add a term +X)i e (^j) 2 to the Hamiltonian, 
which corresponds to a free propagator (X/e)SxY- Finally 



FIG. 6: (Color online) Feynman rules for calculating 1/n cor- 
rections of O(n) model. O(n) indices are omitted. X, Y are 
sublattice indices. The first panel contains free propagators 
and the only vertex in the theory. Straight lines represent 
the spin propagator. The second panel is the one-loop Dyson 
equation Eqn. (|10[) for the A propagator. The third panel is 
the one-loop Dyson equation Eqn. for the spin propaga- 
tor. Thick lines are full propagators. 

we will take the e — > limit. 

The one-loop Dyson equation for the propagator of A 
is shown in Fig. [Sj The inverse of the Green's function 
at one-loop level is 

[G^{k)} XY ={e)S XY - T X , XY = -T X , XY (10) 

where T\ is the self-energy of A (the loop diagram in the 
second panel of Fig. [5]). Here we have taken the e — > 
limit. 

rw-EH) 2 I rf 2 q4°L^(k + q)4°L.yx(q) 

•'BZ 

= -n / d 2 q[M- 1 (k + q)]xr[M- 1 (q)] rx 

There is no summation over sublattice indices X, Y on 
the right-hand side. J BZ d 2 q is the normalized integral 
over the entire BZ. Since the summation over O(n) index 
a becomes a factor of n, the one-loop A propagator is of 
the order 1/n. 

We use this one-loop A propagator to calculate the one- 
loop correction to the spin propagator. 

[Gs?aaQ<)]xY =[M(k)]xy - r S , oa , xy (11) 
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where Ts,xy is the self-energy of spins (the loop diagram 
in the third panel of Fig. 

T S ,aa,XY = H) 2 / d 2 qGf aa XY (k - q)G x , XY (q) 
Jbz 

Again there is no summation over X, Y on the right-hand 
side. 

These integrals cannot be evaluated exactly. Instead 
we use the high-temperature (small (3) expansion to get 
analytical result. We found that up to (i 7 order the one- 
loop correction does not qualitatively change the form of 
Gj^fk). It has similar wavevector dependence of the 
inverse free propagator M(k), therefore the line degen- 
eracy of a > 1 model and the degenerate band of a = 1 
model cannot be lifted at (3 7 order. 

However at /3 s order a qualitatively distinct correction 
appears. The self-energy (the loop diagram) contains a 
term 

/ acos(^f) acos(^) \ 

(l/n)(3 s C acos(^) ^cos(^) 

Vacos(^) ^cos(^f) J 

where fcy = k{ — kj and a constant C = a 2 /(64^/i| c ). 
This looks like a next-neighbor ferromagnetic coupling. 

For a > 1 model we have a degenerate line kx = 
at the saddle point level. This (1/n) (3^ correction will 
favor &2 = 7r which is the M-point in the BZ. For a — 1 
model we have a degenerate band in the saddle point 
approximation. This (l/n)/3 8 correction will favor k\ = 
A)2 = 2n/3 which is the K-point in the BZ, corresponding 
to the -\/3 x v3 spin configuration. 

We notice that a previous high temperature series ex- 
pansion study 3 ^ also lifts the degeneracy of the Kagome 
O(n) model at (3 s order. Their result contains, in some 
sense, corrections to all orders of 1/n, but do not have 
a simple analytical form. Our simpler analytic method 
(expanding in both 1/n and 0) is complementary to their 
linked-cluster series expansion study and our results are 
consistent with theirs in the region of overlap. 

V. QUANTUM LIMIT: EXACT 
DIAGONALIZATION AND SLAVE PARTICLE 
MEAN FIELD THEORIES 

We have performed exact diagonalization on small lat- 
tices of spin-1/2 moments, to study the effect of distor- 
tions in the Kagome model. We used the open source 
ALPS library and applications 31 on an office computer. 
Two different kinds of results are presented; the nature of 
the ground state and low lying excitations, and the ther- 
modynamics (specific heat and magnetic susceptibility). 
The latter requires knowing all eigenvalues of the Hamil- 
tonian, and is hence restricted to small system sizes of 
12 sites (2x2 unit cells) with periodic boundary condi- 
tions. Based on previous studies^ we believe that this 



small system can still produce qualitatively correct high 
temperature properties. For the former, we study system 
sizes upto 24 sites (4 x 2). 

Before discussing the results of exact diagonalization 
for the low lying eigenstates, let us briefly recall the ex- 
pectation from the semiclassical picture developed so far. 

1. For a < 1/2 a colinear ferrimagnetic state with a 
magnetization of 1/2 per unit cell is expected. 

2. For a > 1, the ferrimagnetic chirality stripe state 
is expected, which implies a net spin in the ground 
state and low energy spin excitations at the M 
point as shown in the inset of Fig. [7) 

3. For 1/2 < a < l,wedo not have a firm expectation 
from semiclassics, however, the large-n saddle point 
solution seems to favor a q = state, which would 
also be ferrimagnetic. 

Numerically, we find that the first prediction is remark- 
ably well obeyed even in this extreme quantum limit. On 
decreasing a, the ground state is found to have non-zero 
total spin. Moreover, this is found to happen precisely 
below a = 1/2. The ground state moment is also exactly 
what is expected, for example it is S = 2 for the 2x2 
lattice and S = 3 for the 3x2 lattice. For a > 1/2 the 
comparison is less clear. For example, the ground state 
is a spin singlet on lattice sizes up to 24 sites. However, 
there is a clear tendency of the S = 2 state at the T point 
to drop in energy on moving away from the isotropic 
Kagome point as seen in Fig. indicating perhaps a 
tendency to develop a net moment. On the other hand, 
while the spin gap may be expected to be soft along the 
M point (the wavevector for the chirality state) for a > 1 , 
it turns out that the M point is actually not the loca- 
tion of the lowest spin carrying excitation - which instead 
occurs at different wavevectors; the M' (and equivalent 
M") locations in the case of 2 x 2 system. Similarly, for 
the 4x2 system, the S — 1 excitation energy at the M 
point is higher than those at the M' and M" points (the 
latter two are inequivalent on this lattice). Moreover the 
lowest S = 1 excitation occurs at the M" point. 

It should be noted though that the S = 1 excitation 
energy at the M point decreases rapidly from the 12 site 
to the 24 site latticed and might end up being the lowest 
spin excitation at larger system sizes. Paradoxically, in 
the 1/2 < a < 1 limit, the M point is the location of 
the lowest spin carrying excitation, both in the 12 and 
24 sized systems we studied. We therefore have to leave 
open the question of the validity of the semiclassical 'chi- 
rality stripe' picture in the extreme quantum limit, to 
future systematic numerical studies on bigger systems. 
Finally, we note that as in the isotropic kagome case, we 
have observed singlet excitations inside the spin-gap (the 
energy of the lowest excitation with non-zero spin). 

Thermodynamics: The magnetic dc-susceptibility and 
specific heat results for several different a are presented 
in Fig. [HI In both figures the temperature has been 
rescaled by the average coupling J aV crago = (2 + a)/3 for 
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FIG. 7: (Color online) spin excitation energies at different 
wavevectors from the exact diagonalization study (12 sites: 
2x2 unit cell system with periodic boundary condition). En- 
ergy is measured from the lowest singlet state (ground state 
for a > 1/2). Other spin excitations, S > 2 at any wavevec- 
tor or S — 2 away from the F-point, have much high energies 
than the ones plotted. 
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FIG. 8: (Color online) spin excitation energies at the high 
symmetry points (see inset in Fig.[7]for definition) vs. number 
of sites, for a = 0.75 and 1.25, for 12 site and 24 site systems. 
Due to the geometry of the 24-sites (4x2 unit cell) lattice, 
the M' and M" points are not equivalent. Lines are guides for 
the eye. 



each curve and x is a l so rescaled accordingly. For high- 
temperature (T > 0.2J average ) the dc-susceptibilities for 
different a converge to the a = 1 result. The posi- 
tions of the broad maxima in the specific heat curves 
are also more or less the same for different a. Therefore 
we conclude that the anisotropy does not induce qualita- 
tive difference in these two macroscopic observables for 
high enough temperature (e.g. T > 0.2 J averag0 ). 

Slave Particle Approaches: Other theoretical ap- 
proaches can also be used to attack the problem directly 
from the quantum limit. These methods have been ap- 
plied to the isotropic Kagome lattice and can be utilized 




FIG. 9: (Color online) susceptibility \ an d specific heat Cv 
from the exact diagonalization study (2x2 unit cell system, 
12 spins, with periodic boundary condition). Temperature 
is rescaled by the average coupling J aver agc = (2 + q)/3 for 
each curve. The susceptibilities of different a > 0.5 converge 
to the a — 1 result even at moderate temperatures. The 
positions of specific heat maxima at around T/ Javeragc =2/3 
are consistent between different a values. 



to study the effect of distortion. The Schwinger boson 
technique (large-N Sp(N) approach) has been used to 
study the Volborthite lattice recently 3 - 3 ., where for not 
too large spatial anisotropy the -\/3 x -\/3 state was found 
to persist, although the ordering wavevector is shifted to 
an incommensurate value (the staggered chirality pattern 
remains the same). Fermionic slave particle representa- 
tion of the spins 3 - 4 , as well as the dual vortex formulation 35 
have recently been used to study the isotropic Kagome 
lattice in connection to Herbertsmithite. Extending these 
studies to the Volborthite lattice should be interesting. 
For example, the Dirac fermions in the proposal of Ref. 34 
would remain massless on the distorted lattice as well, 
since the mass term is prohibited by the translational 
and time reversal symmetries that remain intact. 



VI. CONCLUSIONS 

We have studied the distorted Kagome model by sev- 
eral approaches. First we proved that the classical degen- 
eracy is reduced from an extensive one (of the isotropic 
Kagome model) to a sub-extensive one. As a result, 
we found that the ground state ensemble is much less 
connected in the distorted Kagome model compared to 
the isotropic case. One has to change an infinite num- 
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ber of spins (in the thermodynamic limit) in order to 
move from one classical coplanar ground state to an- 
other. This could result in a greater tendency towards 
glassy behavior and may be consistent with the fact that 
spin freezing was observed (not observed) in Volborthitc 
(Herbertsmithite). We then studied the properties of the 
ground state ensemble by enumeration and transfer ma- 
trix methods. Using transfer matrix method we calcu- 
lated the probability of different local spin configurations 
and showed that this consideration may provide an expla- 
nation of the low temperature NMR data in Volborthitc. 

We then studied how this remaining degeneracy can 
be lifted by two novel refinements of various approaches 
to the classical problem. In particular, we used a low 
temperature classical spin-wave expansion to compute 
the effective chirality interactions which lead to a pre- 
ferred ordering pattern. We also studied the large-n O(n) 
model in the saddle point approximation and with 1/n 
corrections, the latter performed in conjunction with a 
high temperature expansion. Our results for the isotropic 
case a = 1 are consistent with previous order-by-disorder 
studies for the isotropic Kagome model, i.e. \/3 x V3 
state is selected. However for a > 1, both classical 
approaches we pursued point to a possible long-range- 
order pattern different from that of the isotropic Kagome 
model. The resulting 'chirality stripe state' doubles the 
magnetic unit cell, has a Fourier component at the M- 
point in the Brillouin zone, and has a net magnetic mo- 
ment (Fig. [1]). Of course, this classical 2D system cannot 
develop a long-range-order at any finite temperature, but 
in the presence of weak inter-layer couplings, the ordering 
pattern we propose is the most reasonable candidate if 
magnetic long-range-order sets in. Exact diagonalization 
studies of small systems showed that the specific heat 
and susceptibility for different values of a do not vary 
much at intermediate temperatures upon the change of 
the anisotropy parameter a. 

Shortly after completion of this work there appeared 
another paper— studying the same lattice but via the 
Sp(N) large-N treatment and perturbation theory. Their 
analysis of the degree of classical degeneracy is in agree- 
ment with our result. 
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FIG. 10: (Color online) brickwall lattice for the transfer ma- 
trix study in Appendix fS] Seven possible single 'brick' con- 
figurations and their probabilities in the thermodynamic limit 
are presented. 



APPENDIX A: TRANSFER MATRIX SOLUTION 

OF THE CLASSICAL GROUND STATE 
DEGENERACY OF THE DISTORTED KAGOME 
MODEL 

In this appendix we derive the asymptotic formula 
of the classical ground state degeneracy in the dis- 
torted Kagome model, and also establish rigorous upper 
and lower bounds to show that the degeneracy is sub- 
extensive. We also study the probability of various local 
hexagon configurations in the ground state ensemble of 
the distorted Kagome model, which is related to NMR 
studies of the VolborthiteiS. 

We stretch the honeycomb chirality lattice horizon- 
tally to make a topologically equivalent 'brickwall' lattice 
(Fig. [IT))) . Chiralities are Ising variables on the vertices. 
For simplicity of derivation we use a different, less sym- 
metric, geometry other than the geometry used for enu- 
meration study in the main text. The lattice consists of 
M rows of 'bricks', each row contains L 'bricks'. We will 
establish the upper and lower bounds, 4 M+i and 2 M+1 , 
for open boundary condition, and the asymptotic formula 
2 M+L for periodic boundary condition in the thermody- 
namic limit. 

It is better to represent the states of the Ising chirality 
variables by domain wall configurations (Fig. [TU)) . As 
in all Ising systems, the number of Ising configurations 
is two times the number of domain wall configurations. 
All possible domain wall configurations within a 'brick' is 
given in Fig. 1101 Number below each 'brick' configuration 
is the probability of that local configuration in a lattice 
with periodic boundary condition in the thermodynamic 
limit, to be derived later. 

There are several important observations: 
(a) If there is a horizontal domain wall crossing one of the 
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vertical edges of the 'brickwall', this domain wall must 
extend through the entire lattice, cutting the entire row 
of 'bricks'. 

(b) The number of vertical domain walls cutting a hor- 
izontal line in the 'brickwall' is conserved from line to 
line. 

(c) Whether there is a horizontal extended domain wall 
in the row of 'bricks' or not completely determines the 
propagation of vertical domain walls from the upper line 
to the lower line. 

(d) If there are two vertical domain walls in the same 
'brick' in the upper line (we call this a 'collision' of two 
vertical domain walls), then there must be a horizontal 
extended domain wall in the row of 'bricks', and we have 
only one choice for the vertical domain wall configuration 
on the lower line. Otherwise for a given vertical domain 
wall configuration on the upper line we have two choices 
on the lower line. 

(e) Vertical domain walls do not cross each other. 

We can obtain an upper bound for the number of chi- 
rality configurations by the following considerations for 
a lattice with open boundary condition, 
(i) The vertical domain wall configurations on the top- 
most line give 2 2L choices; (ii) The horizontal extended 
domain walls give a factor of at most 2 M ; (iii) On each 
row of 'bricks' except for the first row there could be one 
additional Ising degree of freedom depending on whether 
there is a vertical domain wall entering from the top-right 
edge of the rightmost 'brick' (an example of entering ver- 
tical domain wall is given in Fig.[TO]- the third row from 
top). This is at most a factor of 2 M_1 . Combining all 
these factors we get an upper bound 4 M+L for chirality 
configurations on the L x M open boundary lattice. 

We can easily get a sub-extensive lower bound for open 
boundary condition by considering the case that there is 
no vertical domain wall. Then we have 2 M domain wall 
configurations via the M possible horizontal extended 
domain walls. Thus a lower bound of the number of 
chirality configurations is 2 M+ . 

Now we impose the periodic boundary condition on an 
L x M 'brickwall'. Strictly speaking the periodic bound- 
ary condition will introduce two additional non-local con- 
straints on the chirality variables. And it will impose 
constraints on the total number of vertical domain walls 
(must be even) and also horizontal domain walls. They 
are not supposed to change the asymptotic behavior and 
we ignore them for simplicity. 

Define the transfer matrix T xyi where x, y label the ver- 
tical domain wall configurations on the upper and lower 
line of a row of 'bricks', respectively. T xy is the num- 
ber of ways that vertical domain walls in x can propa- 
gate downward to y. Some examples: (i) x is the con- 
figuration where there is no vertical domain wall in a 
line, then the only y satisfying T xy ^ is y — x — 
(no vertical domain wall), and T xx = 2 because there 
could be one, or no, extended horizontal domain wall in 
between, which should be counted as two different ways 
of propagation; (ii) x is the configuration where there are 



vertical domain walls on every edge of the upper horizon- 
tal line, then T xx = 1 because there must be one extended 
horizontal domain wall in between, and T xy = 0, Vy =/= x. 

The number of domain wall configurations is the trace 
of the Af-th power of the 2 2i x 2 2L transfer matrix T, 
which equals the sum of the M-th powers of all eigen- 
values A of T, Tr(T M ) = Ea aM - From the previous 
observations (c) and (d) we have J2 y T xy < 2. Therefore 
all eigenvalues have absolute values smaller than or equal 
to 2. This provides an upper bound 2 2L+AI for domain 
wall configurations. 

Take the thermodynamic limit M — » oo, with L large 
but finite, then the trace Tr(T M ) reduces to the sum of 
the M-th power of the largest eigenvalues (it is 2 and can 
be degenerate), X)a=2 Now we want to construct all 
eigenvectors corresponding to eigenvalue 2. The prop- 
erty of the ground state ensemble is dominated by these 
eigenvectors in the thermodynamic limit. 

Suppose a x is a (left) eigenvector with the eigenvalue 
2j 12 x a xT xy — 2a y . Then we have the following two 
properties: (i) a x > 0, Vx, this comes from the fact that 
T xy > 0; (ii) a x = for x containing a 'collision', this 
comes from the observation (d). 

If there is no 'collision' in x, but there is one vertical 
domain wall crossing the top-left horizontal edge of one 
'brick', and one of its neighboring vertical domain wall 
crosses the top-right horizontal edge of another 'brick', 
we can always bring those two vertical domain walls to- 
gether to make a 'collision', by propagating them down- 
ward (an example is shown in Fig. I10p . Therefore we 
must have a x = for this kind of x, which contains 
both vertical domain walls crossing top-left and top-right 
edges of some 'bricks'. 

Now we can construct all eigenvectors with the largest 
eigenvalue(2) . Take an x containing vertical domain walls 
crossing only the top-left edges of some 'bricks'. It can 
propagate to the next line without change, or shifted by 
one half of the lattice constant. By translating this x 
on the line (with periodic boundary condition) by multi- 
ples of half lattice constant, we find a connected sub- 
space of the state space, denoted by span(a;). Then 
a y = 1, Vy G span(x) is the (not normalized) eigen- 
vector with the largest eigenvalue(2) in this subspace (by 
Perron-Frobenius theorem this eigenvector is unique in 
this subspace). 

The degeneracy of the largest eigenvalues (2) equals to 
the number of distinct subspaces constructed as in the 
previous paragraph, or the number of inequivalent x with 
only the top-left-edge vertical domain walls (inequivalent 
under translation). This is still a non-trivial combinato- 
rial problem, but we have a rough upper bound 2 L and 
a lower bound 2 L /L. Combining all previous consider- 
ations we have the asymptotic form of the number of 
configurations 2 M+L . 

Now we have, in principle, all the eigenvectors relevant 
in the thermodynamic limit. We can find the probabil- 
ities of every 'brick' configuration, or the configuration 
of the six spins in a hexagon in the original distorted 



14 



Kagome lattice. This is related to the 51 V NMR study 
in Bert et al*&, because different local spin configurations 
will produce different magnetic field on the V site. How- 
ever the authors of that experimental paper did not take 
into account the constraints on chirality variables, thus 
their theoretical estimates of the probabilities of different 
local configurations are incorrect. 

First we consider the 'brick' configuration containing 
a 'collision' of vertical domain walls. This corresponds 
to the local \/3 x y3 configuration, which produces the 
largest magnetic field (3 times of a single Cu if a ~ 1, 



in general the factor is 2 + 2a 



on the V site. 



However since our eigenvectors do not contain 'collision' 
the probability of this local configuration is zero. 

Next we consider the configuration where there is one 
vertical domain wall and also one horizontal domain wall 
through the 'brick'. This will produce a smaller magnetic 
field (-s/3 times of a single Cu if a ~ 1, in general the 
factor is \/ (5a — 2)/a 3 ). 

Notice that we have a particle-hole like symmetry. For 
a subspace span(x) discussed in the previous paragraphs, 
where x contains vertical domain walls through some of 
the top- left edges of 'bricks', we can construct another 
subspace span(cc) from a 'complementary' configuration 
x 1 in which there is one vertical domain wall through a 
top-left edge of a 'brick' if and only if there is no vertical 
domain wall through that edge in x. 

Therefore the probability that there is one vertical do- 
main wall through the 'brick' is one half. The probability 
of a horizontal domain wall through the 'brick' is clearly 
also one half for the eigenvectors we consider. Combin- 
ing these two factors we have the probability 25% for 
this type of local configuration. Note that whether the 
vertical domain wall is on the left- or right-side will give 
another factor of one half, hence the 12.5% probabilities 
in Fig. [10] for the two configurations of this type. 

Probability of other configurations can be derived in 
the similar fashion. But all the other local configurations 
will produce very small magnetic field on the V site (for 
a ~ 1). In particular, the configuration with no domain 
wall through the 'brick' has a magnetic field |2 — 2a" 1 1 
times a single Cu field, with the probability 25%. The 
configuration with no vertical domain wall but a hori- 
zontal domain wall has the same magnetic field factor 
|2 — 2a" 1 !, with the probability 25%. The two configu- 
rations with one vertical domain wall but no horizontal 
domain wall have the magnetic field factor \a — l|/a 2 , 
and the total probability 25% (12.5% each). 

Based on these analyses we argue that the 20% slow 
component observed in NMRiS is not due to the local 
•\/3x -\/3 configuration, but rather the configurations pro- 
ducing a smaller (factor rather than 3) magnetic field 
and with a theoretical probability 25% (with one vertical 
and one horizontal domain wall). 



APPENDIX B: DISPERSION OF QUADRATIC 
QUANTUM SPIN WAVE 

In this appendix we present the quadratic (or the so- 
called 'linear') quantum spin wave dispersion of the dis- 
torted Kagome Heisenberg model. We notice that there 
is still a zero-energy band, and the 'spin wave velocity' 
of the dispersive branch vanishes in one direction in mo- 
mentum space. 

We start from Eqn. ^ and do the Fourier transform 
of the bosonic fields, 

b X ,k = N- 1 ' 2 expHk • (R + r x )]&x,R 

R 

where X = A, B, C labels the three sublattices, N is the 
number of unit cells. R are positions of unit cells, Yx are 
positions of the three basis sites within a unit cell, and k 
is the wavevector. 

The quadratic Hamiltonian is then block-diagonalized 



Ho 



k 



• M (k) • ip k + constant 



where ^ = (6^, b Bk , 6^ k , b A _ k ,b B _ k , b c _ k ),^ M(k) 
is a 6 x 6 hermitian matrix, and the summation is over 
the k-points in the BZ. Here M(k) can be written as 



M(k) 




Here P and Q are both 3x3 matrices as shown below, 
and we use the notation c\ = cos(4f ), c 2 = cos(%-), and 
C3 = cos(^j-) with fej = k ■ ej, k 3 = —ki — k 2 - 



Q 



2a + 1 
2a 



(2a-l)c 3 (2a- l)c 2 
4a 2 ci 
(2a - l)c 2 ci 4a 2 

^ c 3 c 2 
c 3 (2a - l)ci 
V c 2 (2a-l)ci 



We need to further diagonalize M(k) by an SU(3,3) 
Bogoliubov transformation. Namely we need an SU(3,3) 
matrix U such that 




and 



U^M(k)U 



w(k) 
v u;(-kl 



where l 3x3 is the 3x3 identity matrix, w(k) is a 3 x 3 di- 
agonal matrix with three branches of spin dispersions as 
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the diagonal elements, because of the inversion symmetry 
w(k)=w(-k). 

In isotropic Kagome model P and Q commute and 
can be diagonalized simultaneously which simplifies the 
calculation. But for general a matrices P and Q do not 
commute. 

A simpler way to get the dispersion is to solve the 
eigenvalues of r • M(k). It is fairly simple to prove that 
the six eigenvalues of r • M(k) are ±c*;,(k), i = 1, 2, 3 in- 
dicating three branches^. The characteristic polynomial 



2/2^ + fix ■ The dispersion is the 



of r • M(k) is 
following 



wi = 0, uj 2 ,3 = y/ 2 T Va 
where A = /f — /4 and 
/ 2 =2a 2 + 1 - 2a" 1 + 2a" 2 - (2a 2 - 1) cos(fci) 
— a" 1 [cos(fc2) + cos(/c 3 )] 
,(«-l) 2 



A =2- 



x {2 + a 2 + a 2 cos(fci) - 2a[cos(fc 2 ) + cos(/c 3 )]} 



Although the dispersion has become much more compli- 
cated than the Kagome case, the zero-energy band still 
exists. 

When a = 1, /f — f± = 0, we have w 2 = oj 3 = 
y3 — cos(fci) — cos(fc 2 ) — cos(/c 3 ). For small |k| the dis- 
persion becomes w 2 = w 3 ~ v^i + ^i^ 2 + fc 2 oc |k|. 
Thus we have two 'linear' spin wave branches. 

However, as long as a / 1, we have o> 2 7^ w 3 and 
cj 3 (k = 0) = 2 1 1 - a _1 | > 0. We still have one Goldstone 
mode because w 2 (k = 0) = 0. But the small wavevector 
dispersion is drastically changed, w 2 ~ \J (a 2 — l/4)fc 2 cx 
|fci|. Namely the 'spin wave velocity' in the fc 2 direction 
(vertical direction in k-space) vanishes. 



APPENDIX C: CLASSICAL SPIN WAVE: 
QUADRATIC THEORY AND CHIRALITY 
INTERACTIONS 

Let us start from Eqn. ([5|) , replace e y and e z by e v and 
e z , and do the Fourier transforms of e v and e 2 (see the 
previous appendix for notation) 



%,k = cx PH k • (R + rx)]^ R 

R 

%,k = N- 1 / 2 cx PH k • ( R + rxWxv 

R 

The quadratic Hamiltonian can be block-diagonalized 

^! = E^(k)x k 

k 

#1 = E4^(k)<£k 




FIG. 11: (Color online) chirality-chirality couplings calculated 
here. Equivalent couplings under space group symmetry are 
not shown. 



TABLE II: Effective chirality couplings Ji... 6 (see Fig. [TTj) 
divided by T 2 , for J z = 0.1, J y = 0.01. Positive number 
means antiferromagnetic coupling. 



a 


Jlv 


Jlh 








1 


0.9702 


0.9702 


0.2614 


0.2614 




1.5 


1.8231 


-0.3340 


-0.01294 


0.8895 




a 






J 4u 


J 4 m 


J id 


1 


0.1916 


0.1916 


0.002661 


0.002661 


0.002661 


1.5 


0.4318 


0.3631 


0.04897 


-0.2770 


0.007420 


a 






J (i v 


Jfih 




1 


0.002924 


0.002924 


0.002914 


0.002914 




1.5 


-0.06089 


0.01172 


-0.03724 


0.1996 





where 

(^A,-k> ^B,-k' ^c, 



(^A.-k^B.-k^a-k) and 4 = 
_ k ), and Mj, iZ (k) are both 3x3 
matrices, shown below, where we use the notation 
ci = cos(^ L ), c 2 = cos(^-), and c 3 = cos(^). 




My(k) = 



M,(k) 



It is straightforward to check that M z (k) has a 
zero eigenvalue with (not normalized) eigenvector 
(asin(fci/2),sin(A: 2 /2),sin(fc 3 /2)) for all k; and M y (k = 
0) has a zero eigenvalue with eigenvector (1, 1, 1). 

For small |k| the dispersion of the lowest branch of 
M y (k) is approximately (l/6a)(a 2 fc 2 + fci/c 2 + /cf). 

Now we consider the calculation of the chirality in- 
teractions. Each chirality interaction is calculated by 
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thirty-six terms in (H3) 2 , we show here an example in 
Fig. O Chiralities 771 and 775 arc defined on triangles 
ABC and DEF in the distorted Kagome lattice, respec- 
tively. 771 determines the sign of the angles between spins 
on ABC sites, 9 AB = 7/16*0, BC = -2t7i#o, &ca = Vi^o- 
6*de, #ef, Ofd arc determined in the similar way by 775, 
and 9ji — —6ij. Plug these into Eqn. ([5]), then the rele- 
vant terms in (H 3 ) 2 are 2i]iri 5 (h AB + h BC + /jca)(/jde + 
/ief + /ifd) where 

h AB = sin(6» Q )(e^ B - e B 60 
/iBC=sin(-2e )(^c-e^B) 
h C A = sin(6>o)(e^£ A - e^c) 

Here we use & = [(ef) 2 + (ef) 2 ], — arccos(— l/2a) 
and /iDE, ZiEFj ^fd are obtained by replacing subscripts 
ABC by DEF respectively. 

According to Eqn. © the ef- 
fective chirality-chirality coupling is 



-T 2 ((h AB + h BC + h C A)(hr>E + h E F + %d)) - Ex- 
panding this expression, we have thirty-six terms, each 
of the form {^((.j^j^m) which can be further expanded 
into four terms (e^)^^) 2 ^ + (^(^) 2 ?I(^) 2 ) + 

m^k^rnf) + mi??k(Z Z m )\ Each term in the 
last expression can be expanded into a sum of prod- 
ucts of three two-point correlators by Wick theorem. 
The two-point correlators are computed following the 
standard routine in all quadratic theory, e.g. 

(^,o^,r>o = / rf 2 k[M- 1 (k)] AB e ik -( R + r B- A ) 

for the A-sublattice site in the unit cell at origin and 
the B-sublattice site in the unit cell at position R. We 
calculated up to the sixth neighbor chirality couplings 
(Fig. QT]) . Some data are presented in Table [Til 

All the above mentioned calculations in Appendix [B] 
and Appendix[C]were done by the software Mathcmatica. 
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